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J> ■ Abstract 

This paper is an attempt to present and discuss at some length the Singular Mani- 
oc ■ fold Method. This Method is based upon the Painleve Property systematically used 
as a tool for obtaining clear cut answers to almost all the questions related with 
Nonlinear Partial Differential Equations: Lax pairs, Miura, Backlund or Darboux 
| Transformations as well as r-functions, in a unified way. Besides to present the ba- 
sics of the Method we exemplify this approach by applying it to four equations in 
(1 + l)-dimensions. Two of them are related with the other two through Miura trans- 
formations that are also derived by using the Singular Manifold Method. 

43 

1 Introduction 



\ 1.1 Integrability and the Painleve Property 

X 

The beginnings of the study of singularities in the complex plane for differential equations 
has always been attributed to Cauchy j^]. Cauchy's main idea was to consider local 
solutions on the complex plane and to use methods of analytical prolongation to obtain 
global solutions. For this procedure to work, a complete knowledge of the singularities of 
the equation and its location in the complex plane is required. In this sense, it is essential 
to distinguish between two types of singularity. 

• Fixed singularities: Singularities determined by the coefficients of the equation 
and its location does not therefore depend on initial conditions. 



Movable singularities: Singularities whose location on the complex plane does 
indeed depend on the initial conditions. 



Copyright ©1998 by P.G. Estevez, E. 
Conde and P.R. Gordoa 



Miura, Backhand and Darboux Transformations 



83 



The truly decisive step towards the elucidation of the relationship between the analyti- 
cal structure of a system and its integrability is attributed to the Russian mathematician 
Sofia Kovalevskaya ^4|. Her work focused on the study of the motion of a rigid solid 
with a fixed point from an analysis of the singularities of the solutions ^4|. Kovalevskaya's 
work was completely new and also addressed to uniquely determine the parameter values 
for which the only movable singularities of the solutions on the complex plane were poles. 

Although Kovalevskaya's work was apparently afflicted by a lack of followers, in the 
last decade of the nineteenth century some mathematicians focused their attention on 
the classification of ordinary differential equations (ODEs) on the basis of the type of 
singularity their solutions were able to exhibit. 



It was the French mathematician Paul Painleve [32] who, following the ideas of Fuchs, 
Kovalevskaya, Picard and others, completely classified first order equations and studied 
second order equations. In this last case, he found 50 types of second order equations 
whose only movable singularities were ordinary poles. This special analytical property 
now carries his name and in what follows will be referred to as the Painleve Property 
(PP) |35|. Of these 50 types of equation, 44 can be integrated in terms of known functions 
(Riccati equations, Elliptic Functions, Linear Equations, etc) and the other six, in spite of 
having meromorphic solutions, do not have algebraic integrals that allows one to reduce 
the equation to quadratures. Today these are known as Painleve Transcendents. The 
main contribution of Paul Painleve lies in that he established the basis for a theory that, 
unlike what had been believed until then, would allow one a priori, by singularity analysis, 
to decide on the integrability of an equation without previously solving it. 

Although there is no definitive proof of why singularity analysis for an equation turns 
out to be a test of integrability, some circumstances do seem to corroborate this. For 
example, it has been found that Painleve Transcendents often appear in similarity reduc- 
tions of equations with solitons Q , |2j| . It is also intriguing to note that always a certain 
relationship seems to exist between equations with the PP and Isomonodromy Transfor- 



mations of certain linear equations [14]. The validity of Painleve's analysis as a suitable 
procedure for detecting integrability could be related to the combined study of Algebraic 
Geometry and Abelian Function theory |j4|. Jacobi elliptic functions have associated 
with them a genus 2 Riemann surface (torus with a hole) and on this surface all algebraic 
curves are meromorphic functions. Although hyperelliptic functions cannot generally be 
parametrized in terms of meromorphic functions, Jacobi observed that certain combina- 
tions of hyperelliptic integrals do have meromorphic inverses (this, for example, is the case 
of the integrals obtained by Kovalevskaya in the fourth integrable case for the rigid solid) 
and these are called Abelian integrals. Simultaneous study of these integrals together 
with the associated Riemann surfaces could be crucial for establishing the final proof of 
a reliable test of integrability based upon the ideas first developed by Paul Painleve and 
Sofia Kovalevskaya. 

1.2 The Ablowitz, Ramani and Segur algorithm 

Between 1955 and 1960 the Korteweg-de Vries equation reappeared in the work of Fermi- 
Ulam-Pasta Pa and in the context of plasma physics. Towards 1965, with Zabusky 



and Kruskal [39], the concept of soliton emerged for the first time. The Soliton was 



an entity describing solitary wave solutions interacting among themselves without any 
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change in shape except for a small change in its phase. With this discovery in mind the 
Inverse Scattering Technique (1ST) was developed ||, initially allowing one to solve the 
KdV equation and then many integrable equations with soliton solutions. The incredible 
success obtained with the 1ST technique prompted to Ablowitz, Ramani and Segur Q to 
develop an algorithm (similar to that used by Kovalevskaya in the rigid solid problem) to 
determine whether an ordinary differential equation had the PP. As mentioned above, an 
ODE is said to have the PP if the only movable singularities of its solutions are poles. This 
same property can be stated by saying that all solutions are singlevalued except in the 
fixed singularities of the coefficients. The ARS (Ablowitz, Ramani and Segur) algorithm 
is a method for determining the nature of the singularities of the solutions of an ODE on 
the basis of an analysis of their local properties. 

Until now we have considered the study of singularities within the context of systems 
described mathematically by ODEs. In view of the results obtained with the ARS al- 
gorithm, it seems natural to ask ourselves whether some other method for dealing with 
differential equations in partial derivatives (PDEs) could in principle be available. In 
particular, to find a new version of the PP that can be applied directly to the study of 
integrability for PDEs could indeed be extremely useful. One important problem here is 
that the solutions of a PDE are functions of at least two independent variables and their 
analytical continuation is clearly more complicated than in the case of ordinary equations. 

The study of similarity reductions for PDEs that can be solved by 1ST led Ablowitz, 
Ramani and Segur to formulate what has now become know as the ARS conjecture: 
"Every ordinary differential equation that can be obtained as the similarity reduction of 
a PDE solvable by 1ST has the PP up to a smooth change of variables". This conjecture 
provides a necessary condition for checking whether a PDE is integrable or not. Ablowitz, 
Ramani and Segur [|| and Mc Leod and Olver [^9| have tested some weak versions of this 
conjecture. Such demonstrations are based on the fact that if a PDE can be completely 
integrated its solutions can be evaluated in terms of the Gel'fand-Levitan-Marchenko linear 
integral equation that appears in the 1ST. 

The ARS conjecture can therefore be understood in the sense that if it is possible to 
reduce a PDE to an ODE that does not have the PP (even after a suitable transformation 
of variables) it may be concluded that the PDE is not integrable. An algorithmic procedure 
has recently been put forward for determining similarity reductions for PDEs. The essence 
of the procedure is the study of the Lie symmetries. To check that a PDE has the PP 
using the ARS conjecture one must find all the possible similarity reductions and check 
that all the resulting ODEs do have the PP even though one has to make transformations 
of variables. In this context the ARS conjecture is on the one hand tedious (owing to the 
huge number of reductions to ODEs shown by some equations) and on the other hand 
becomes less useful as the number of symmetries shown by the equation decreases. It is 
also not clear which transformations of variables are permitted when checking whether the 
corresponding ODE is of the Painleve type. In particular for equations that do not have 
symmetries the ARS conjecture is quite useless as it is not possible to obtain similarity 
reductions from usual group-theory procedures. 

The obvious limitations of this method suggest that it would be interesting to have 
available a direct method (analogue of the ARS algorithm for ODEs) that would allow 
one to decide whether the PDEs under study are integrable. In this way one is extending 
the definition of the PP to partial derivatives using the original idea of Painleve and 
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developing an algorithmic method to determine whether the equations have this property 
(without the requirement of considering all their similarity reductions to ODEs) and hence 
to decide whether they are integrable or not. 

1.3 The Weiss, Tabor and Carnevale algorithm 

The main difference between analytical functions of one and several variables is that the 
singularities of the latter are not isolated. If f(z\, . . . , z n ) is an analytical function of n 
complex variables z% (i = 1, . . . , n), the singularities of / are in manifolds of (2n — 2)- 
dimensions. These manifolds are determined by conditions of the form 



where x 1S an analytical function in a neighborhood of the manifold defined by (1.1). 
When this manifold depends on the initial conditions it is called a movable singularity 
manifold. The existence of these singularity manifolds suggests the need for introducing 
the PP concept for PDEs. This leads to a procedure to check whether the equations have 
such a property in a way that makes possible to evaluate the question of integrability by 
a unified analysis of singularities for both ODEs and PDEs. This was the work carried 
out by Weiss, Tabor and Carnevale (WTC) |35f| . According to these authors, we say that 
a PDE has the Painleve property (PP) if its solutions are singlevalued in a 
neighborhood of the manifold of movable singularities. 

The WTC method also allows one to successfully apply some of the techniques de- 
veloped for integrable systems to systems that are not completely integrable. Certain 
restrictions may be imposed on x or on the parameters of the equation such that the 
solutions thus obtained only have poles as movable singularities. In this case, the singu- 
larity manifold x is no longer an arbitrary function and the equation is said to have the 
conditional PP Pj and hence is partially integrable. 

1.4 The singular manifold method 

Weiss fl36| , |37[ introduced the singular manifold method (SMM) which is an efficient 
algorithmic method to find the typical properties of integrable systems. If a PDE has the 
PP we have seen that its solutions can be expressed as a Laurent series in the form 



It is possible in any case to truncate the expansion series at a certain term in order to 
obtain particular solutions of the equation. If the expansion is truncated at the constant 
term (understood as the one that goes with cjp where we use 4> instead of x f° r the manifold 
in the truncated expansion), expression (1.2) reduces to: 





oo 




(1.2) 



3=0 



u = u (j) a + ui(j) 1 a H Vu a . 



(1.3) 



It is interesting to note that some attempts have been made in order to truncate the 
Laurent series at higher orders [|33[|; however here we shall only consider truncation at the 
constant term. Substitution of (1.3) in the corresponding PDE leads to an overdetermined 
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system of equations for cf), Uj and their derivatives. The essential point is that the sin- 
gularity manifold <fi is no longer an arbitrary function but rather -as we now shall see- it 
must fulfill certain equations due to the truncation condition. That's why the truncation 
of the Painleve series is the basis of a method called Singular Manifold Method (SMM) 
that has been proved to be extremely successful in studying nonlinear PDEs. Many of the 
properties of such equations can be obtained through the SMM. Let us summarize some 
of them. 



The truncation (1.3) of the Painleve series has itself the meaning of an auto-Backlund 



transformation between two solutions of a PDE [37], [12] 



The Lax pair can be obtained through the Singular Manifold equations pq ], [30]. 
The relation between the singular manifold method and nonclassical Lie symmetries 



of the truncated solutions has been studied in [10]. 



The method of Hirota g(J, @, @ is known ciS db powerful procedure for generating 
multisoliton solution for PDEs. It essentially consists in bilinearizing the differential 
equation by an ansatz reminiscent of the Painleve truncated expansion. 



The WTC method also provides an iterative procedure for generating solutions [17] 
from the Lax pair and from the corresponding auto-Backlund transformation, where 
the corresponding singularity manifold <p is determined in each step and after n steps 
the solution can be expressed in terms of the product 0i, (f>2, ■ ■ ., <p n from which it 
is then possible to construct the Hirota r function associated with the solution with 
n solitons. The relationship between singular manifold and Hirota's r-functions ||, 
[13 1, [ll(| has also been clearly established. 



11 



The Darboux transformations of a PDE [pq] , wm, ||, [27] are also an important 
procedure to obtain solutions of PDEs. The connection between SMM and Darboux 



transformations has been explained in different references [11 1 



1.5 Plan of the paper 

After the previous glimpses of evidence in regard to the relationship between the PP and 
the integrability conditions various directions in the search for integrable PDE become 
evident. On the one hand the Painleve test does identify integrable systems and on the 
other hand the Singular Manifold Method appears as a systematic technique for finding 
Backlund and Darboux transformations, Lax pairs, Soliton Solutions etc. Our point of 
view -as we shall show henceforth- is that the PP should be used not only as a test of 
integrability but also as a fruitful source of information of practically all the important 
features of the Non Linear Partial Differential Equations. This paper tries to go an 
step further in this direction, adding the Miura transformations to the above mentioned 
properties arising from the SMM. 

Sections 2 and 4 deal with the application of the SMM to two equations in (1 + 1)- 
dimensions as the AKNS (Ablowitz-Kaup-Newell-Segur) H and non local Boussinesq 
equation NLBq [|26[], |38| | equation. The Lax pair, Darboux transformations and Soli- 
tonic Solutions are thus fully obtained in these cases. This also shows that our analysis 
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becomes not only a conceptual piece of information but also an algorithmic tool that can 
be systematically used. 

Sections 3 and 5 are devoted to the study of the ShG (sinh-Gordon) || and KS (Kaup 
system) [^] systems. The Miura transformations between these equations and AKNS 
and NLBq respectively are obtained by using the SMM. Backlund transformations for 
AKNS and NLBq are derived. The two component induced Lax pairs for ShG and KS are 
identified by the same procedure. 

Section 6 is one of conclusions. Some lengthy and/or auxiliary calculations are relegated 
to Appendices A to E. 



2 The AKNS equation in (l+l)-dimensions 

The well known 0] AKNS equation in 1+1, 

= M yxxx + 4M y M xx + 8M x M xy (2.1) 

is a nice and easy example to start to show how the method works. 

The PP for this equation means that all solutions of (2.1) can be written as a series of 
the form (see |35|1): 

oo 

M = Y J M,x 1 -\ (2-2) 

where a is the leading index and x is an arbitrary function of x and i, depending on the 
initial data, that is usually called singularity manifold. Mj are analytical functions 
of t and x in the neighborhood of x = 0. It is a trivial exercise to substitute (2.2) into 
(2.1) and to check that if the leading index is a = 1, the series (2.2) satisfies (2.1) for any 
functional form of x- It could be said that (2.1) has the PP. 



2.1 Truncated expansion. Backlund transformations 

The SMM is based upon the above defined PP. It requires the truncation of (2.2) at 
the constant level j = a. It means, for our equation (2.1), that the truncated solutions 
should be: 

M' = M+— , (2.3) 

9 

where we have called M = M\ and M' to the truncated solution. We have used here (ft for 
the singularity manifold, instead of x 5 to emphasize that the truncation implies that the 
singularity manifold is no longer an arbitrary function but a function that is "singularized" 
by the fact that it should satisfy some definite equations as we will see later. We shall 
be calling this manifold (ft singular manifold p6|], p7| ] and the method based on the 
truncation of the Painleve series is the so called Singular Manifold Method. 

The substitution of the truncated expansion (2.3) into the equation (2.1) provides the 
following results (see Appendix A): 

• M as well as M' should be solutions of (2.1). It means that (2.3) could be consided as 
an auto-Backlund transformation between two solutions M' and M of the same equation. 
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The solution M can be written in terms of the singular manifold in the following way: 

Mx = ~(\) (^ + y + 2A )> ( 2 - 4 ) 



M y = -(-v y + 2Xq). (2.5) 
The notation that we have used is [Oil: 



v = -^, (2.6) 

Q = % (2-7) 

and A is an arbitrary constant that, as we will see, plays the role of the spectral parameter. 

• The singular manifold equations. The equations that the truncation procedure 
implies for 4> are: 

s y = 4Ag x , (2.8) 
where s is the schwartzian derivative defined as: 

v 2 

s = v x -—. (2.9) 

Furthermore the compatibility condition <f> xx t = 4>txx between the definitions (2.6) and 

(2.7) requires: 

v y = (q x + qv) x =J> s y = q xxx + 2sq x + qs x . (2.10) 

It is not difficult to prove that the singular manifold equations are nothing but the AKNS 
system once again. In fact with the change of variables: 

s = i Px + 2X, (2.11) 
q = &. (2.12) 

(2.8) is trivially fulfilled and (2.10) is written as: 

= p vxxx + kpyVxx + 8p x p xy (2.13) 
that is obviously the AKNS system. 
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2.2 Lax pairs 

As we have seen above, the singular manifold equations, written in terms of v and w are: 
v xy - vv y = A\q x , (2-14) 

v y = {q x + qv)x (2.15) 

that can be considered as a new system of nonlinear equations. If we apply the Painleve 
analysis to this system, the leading terms are (using ip for the singularity manifold). 

v ~ VQip a , q ~ qoil> h . 

The substitution in (2.14-15) provides: 

a = -1, b = -2, v = 2tp x , q = --i/} x ip y . 



These leading terms provide the key for the linearization of the truncated solutions (2.3-4). 
Actually if we substitute v for its dominant term. 

v = 2% => 4> x = i> 2 - (2.16) 

w 

In such a case (2.4) is: 

= ifrxx + (2M,; + X)ip (2.17) 
and from (2.16), (2.5) and (2.10) we obtain: 

thy 1 

or 



2-f- = q x + qv = —(2M xy + v xy + 2vM y + vv y ) 



= 2\%j) y + M xy i/) - 2M y ip x . (2.18) 

(2.17) and (2.18) are precisely the Lax pair for AKNS. To summarize it is possible to say 
that the Lax pair is nothing but the singular manifold equations in which the 
eigenfuctions are directly obtained from the singular manifold through (2.16). 

2.3 Darboux transformations 

Following an idea of Konopelchenko and Stramp [23], we can consider the Lax pair itself 
as a pair of coupled nonlinear equations between M and V- Let us now explain how to 
proceed. 

As far as M' is also a solution of (2.1) an associated singular manifold (j)' 2 linked to an 
spectral parameter A2 can be defined just by defining 

4>2 X =^2 (2-19) 
a Lax pair for M' can be written as 

= ^ 2xx + (2M^ + X 2 )^ 2 , (2.20) 
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= 2A 2 ^ + M' xy ip 2 ~ 2M^' 2x , (2.21) 

where the notation means that ip' 2 1S an eigenfunction for M' with eigenvalue A 2 . If we 
call <pi and 4> 2 two singular manifolds for M attached to spectral parameters Ai and A 2 
respectively the corresponding eigenfunctions are defined as 

cj> Xx = i>l (2.22) 

<t> 2x = V|, (2.23) 
and the Lax pairs take the form 

= ^l« e + (2M a . + Ai)^i, (2.24) 

= 2Ai^i» + M xy ipx - 2M y ^ lx , (2.25) 

= thxx + (2M x + \ 2 )ih, (2.26) 

= 2\ 2 4> 2y + M xy ip 2 - 2M y i> 2x . (2.27) 
If we use the singular manifold <f>i to construct the truncated Painleve expansion 

M' = M + — (2.28) 
91 

and we then look at (2.20-21) as a system of nonlinear coupled equations a similar expan- 
sion should be performed for tp' 2 . That is to say: 

</4 = ^2 + -^. (2.29) 

The substitution of the truncated expansions (2.28-29) in (2.20-21) provides the func- 
tional form for 0. The result is 

e = -^in(^i,^j), (2.30) 

where 



fifth, fh) = [ Xl -\ 2 J ^ 2x ~ Mlx) ' (2,31) 

The expansions (2.28-29) can be considered as transformations that leave invariant the 
Lax pair (2.18-19). In this sense they are Darboux transformations. It should be pointed 
out that we use the singular manifold <$>\ to actually realize the transformation (2.28) but 
not the eigenfunction tp\ as it is usual [28] in the Darboux transformations. Nevertheless 
eigenfunctions and Singular Manifolds are trivially related through (2.22) and therefore: 
With two eigenfunctions ipi and ip 2 for M, we can construct an eigenfunction 
ip' 2 for the iterated solution M' . That is why we call them Darboux transformations. 
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2.4 Hirota's function 

Furthermore, (2.19) is a nonlinear equation that relates 2 and ip' 2 - ^ means that the 
singular manifold 0' 2 itself could also be expanded in terms of 0i 

2 = 02 + ^ (2.32) 

and by substituting this expansion in (2.19) we obtain: 

A = -[ft(Vi,V2)] 2 . (2.33) 

The procedure described above could be easily iterated. The singular manifold 2 for M' 
can be used to construct a new solution 

M" = M' + ^ (2.34) 

02 

that combined with (2.28) can be written as: 

M" = M + —, (2.35) 

where 

712 = ^201 (2-36) 
and by using (2.32) and (2.33) 

ria = 0201 - [0(^i, V2)] 2 . (2.37) 

It is an interesting point to note that the function ti 2 for the second iteration is not 
a Singular Manifold but it can be constructed from two Singular Manifolds of the first 
iteration. The SMM provides algorithmically, but sometimes considered just a clever 



ansantz [J16|J, [18], the Hirota's bilinear method [19|. It also provides the way to construct 



solutions for the r-function, as we will see in the next section. 
2.5 Solitonic solutions 

The easiest nontrivial solutions can be obtained from the seminal solution 

M = a y. (2.38) 
For this solution, exponential solutions of (2.24-27) are 

ipi = exp ^kiX - -jr.y^j ' ( 2 - 39 ) 

where 

Xi = -k\ (2.40) 
and the corresponding manifolds are 

l = ^-(« J + V' 2 ), (2.41) 
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where on are arbitrary constants. The equation (2.31) provides: 

n^l^H^-j^M* (2-42) 
and (2.37) also yields 

ri2 = TlV(«i + ^i 2 )(«2 + ^l) - J l ft\2 - ( 2 - 43 ) 
4fcifc 2 {ki + k 2 y 

We can write the first and second iteration as: 

M' = a y + ^, (2.44) 

M" = a y + ^, (2.45) 

7"! 2 



where 



f-d + FO, (2-46) 



ri2 = {1 + F x + F 2 + A 12 F 1 F 2 } , (2.47) 



aj = exp (2kiX . 

T* 

fc 2 



2 



exp [2k i [x-^y-x 0i )) t (2.48) 



(2.44) corresponds to the one-soliton solution (see Fig. 1) and (2.45) to the interaction of 
two solitons (see Fig. 2). 

3 The sinh-Gordon equation. 

Miura Transformation to AKNS 

We shall be analyzing in this section how the SMM and the Singular Manifold equations are 
able to provide also information about Miura transformations between nonlinear PDE's. 
In particular, this section is devoted to obtain the Miura map between AKNS and the 
sinh-Gordon equation. 

The sinh-Gordon equation || 

U xy = sinh 2U (3.1) 
can be written as the system 

Q = u xy + 2uT)y, (3.2) 

= Vx + u 2 (3.3) 
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through the change 

u = U x , (3.4) 
rj y = — cosh 2U. (3-5) 



This system has the same problem as the sine-Gordon equation [31|. It has two Painleve 



branches. In fact if we write the solutions of (3.2-3.3) as a Painleve series 



oo 

u 



j=0 

oo 

v = J2vjX J ~ b - 

j=0 

The leading indexes are a = b = 1 but the dominant terms are: 

n = ±Xx, (3.6) 

m = Xx- (3-7) 

The ± sign of uq means that there are two possibilities for the expansion. This may seem 
at first a problem when one attempts to apply the SMM since if we choose a definite sign 
in the expansion we will be loosing information about the equation. This problem has 
been discussed during the last years (see ||, Q|, ||). In these papers a modification 
of the SMM appears to be necessary. We need in fact to deal with two Singular Manifolds 
altogether. Obviously this fact represents a non trivial complication for the calculations. 

We present here the easiest form to work with these two Singular Manifolds. The idea 
is the following: The two branches (3.6-7) suggest the following set of changes for the 
functions: 

u = m — rh, (3-8) 

rj = m + rh (3-9) 

in such a way that m and rh should have only a Painleve branch. It is now necessary 
to look for the equations that m and rh should satisfy. For this purpose we introduce in 
(3.2-3) the change (3.8-9). Adding and subtracting the result we obtain: 

m xy + 2(m — rh)m y = 0, (3.10) 
rh xy — 2{m — rh)rh y = 0. (3-H) 



From these equations we also obtain the additional information that we discuss in the next 
Subsections. 
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3.1 Miura Transformation 

The rh can be obtained from (3.10) and substituted in (3.11). The result is that m obeys: 

= 2m y m XX y + 8m x m 2 - m 2 xy (3.12) 
that is the integrated version of the AKNS equation 

= m yxxx + im y m xx + 8m x m xy . (3.13) 
In a similar form, we can obtain m from (3.11) and the substitution in (3.10) is 

= 2ihyih xxy + 8rh x rh 2 - rh xy (3.14) 
that it is again the integration of the AKNS equation 

= rh yxxx + 4rh y rh xx + 8rh x rh xy . (3.15) 

As a consequence of what has just been said the change of functions (3.8-9) can be inverted 
to yield 

2m x = u x + ri x = u x - u 2 , (3.16) 

2rh x = i] x — u x = —u x — u 2 (3-17) 

which represents the Miura transformations between the sinh-Gordon system (3.2-3) and 
the AKNS equations (3.13) and (3.15). 

3.2 Backlund Transformations 

As an additional result we obtain that the solutions (u, r/) of sinh-Gordon can be con- 
structed by using two solutions m and rh of AKNS. Nevertheless these solutions are not 
independent. They are related by (3.10-11) that can be written as: 

m = m-\ -, m = m -\ (3.18) 

2 m y 2 m y 

(m y rh y ) x = 0. (3.19) 

These equation (3.18-19) can easily be recognized as Backlund transformations between 
the two solutions m and rh of AKNS. 

3.3 Singular manifold method with two manifolds 

From the above discussion it is easy to understand why it seems reasonable to talk about 
two Singular Manifolds, one for the expansion of m and the other for the expansion on 
rh. Let us call <f) the singular manifold for m and (f> the singular manifold for rh. The 
truncated expansions are 

m ' = m + ^, (3.20) 
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m' = m + ^, (3.21) 
and the corresponding expansions for u and rj are: 

/ . 4>X fix rr, n r>\ 

u = u-\ — — , (3.22) 

4> (j) 

V' = ?? + % + ^- (3.23) 

However <fi and <fi are not independent because m and m are related by the Backlund 
transformation (3.18-19). This is reminiscent of the requirement that u and r\ satisfy (3.2- 
3). In fact substituting (3.22-23) in (3.3) we obtain the coupling condition between 4> and 
4> (see Appendix B) 

^k = A^ + lk, (3.24) 

4> 4> 4> 4> 

where 

A = l + u, (3.25) 

A= V --u. (3.26) 
The notation is the one defined in (2.6). 
3.4 Lax pair for sinh-Gordon 

The derivative of (3.24) with respect to x (see Appendix B) provides: 

A x = A{v-A-A) = V -=^, (3.27) 

A x = A{v-A-A) = V -^. (3.28) 

We should remember at this point that (f> and <f> are Singular Manifolds for AKNS and 
that the change 

4>x = *P 2 , k = V> 2 (3.29) 

relates the Singular Manifolds with the eigenfunctions of the Lax pair (2.17-18) of AKNS. 
By combining (3.29) with (3.27-28) we can integrate out the variable x, and finally obtain 

A = a%, (3.30) 

W 

A = at, (3.31) 
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where a and a are constants. (3.25-26) can be now written as: 

ip x = aip — uip, (3.32) 

i\) x = aip + uip, (3.33) 

where tp and ^ are solutions of the Lax pair of AKNS. That is equivalent to 

= il> xx + (2m, + A)V, (3.34) 

= 2\ip y + m xy ip — 2m y ip x , (3.35) 

= 4* + ( 2 ™z + A)^, (3.36) 

= 2\i(jy + m xy ip — 2rh y ip x . (3.37) 
The compatibility between (3.32-33) and (3.34), (3.36) requires 

A = A = -aa. (3.38) 
The y component of the Lax pair is given by (3.35), (3.37) and can be written as 

2aip y = -(u y + r) y )i£>, (3.39) 

2aijjy = (u y - r]y)ip. (3.40) 
Therefore, the Lax pair for sinh-Gordon can be written in its usual matrix form as: 

(3.41) 



t) = ( ~ U a ) ( t 



\ - f u y + v y \ ( t 

2a ^ /„ ~ \ u y-Vy J \ $ 



(3.42) 



Conclusion: The method that we derived in Section 2 to obtain Darboux transformations 
and solutions for AKNS can be directly applied to sinh-Gordon. The construction of 
solution of sinh-Gordon can be done through (3.8-9) by using two solutions of AKNS 
related by the Backlund transformation (3.18). From the point of view of the Singular 
Manifold this implies that 4> and <fi are related by the coupling condition (3.24). Using 
(3.29-31) this condition can be written in a much easier form as: 

a4> + &</) = i/nf>, ( 3 - 43 ) 

A = A = -aa. (3.44) 

We believe that the splitting (3.8-9) is the key to solve the long standing problem that 
concerns to the application of the SMM to equations with two branches. In the next 
sections we will return to the same topic in a different but not unrelated context. 
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4 Non local Boussinesq equation 

The following system of equations 

N x = M t , (4.1) 

M x N t = M X M XXX + 2Ml + M 2 - M 2 XX (4.2) 

or equivalently 

Ml{M tt - M xxxx ) = AMlM xx + 2M x (M t M tx - M XX M XXX ) - M XX (M? - M xx ) (4.3) 

has been considered in pf|, [ 38 1 as related with the Kaup (sometimes called Classical 
Bousinesq) system through a Miura transformation. The Kaup system is a good example 



of system with two branches [21], ||, ||. Nevertheless (4.1-2) has only one branch. That 
is why we are interested in this section in studying the relationship of both equations 
from the point of view of the SMM and to relate the corresponding results with the Kaup 
system in Section 5. 

4.1 Truncated expansion. Backlund transformations 

The leading terms of (4.1-2) are 

M~ XxX -\ 

N ~ XtX' 1 
that suggests the truncated expansion 

M' = M + ^ —>N' = N + ^. (4.4) 

From this equations we obtain the following set of results (see Appendix C): 

• (4.3) could be consided as an auto-Backlund transformation between two solutions 
M' and M of the same equation. 

• The solution M can be written in terms of the Singular Manifold as: 

M *= (i) ((™ + 2A) 2 -t; 2 ), (4.5) 

M t = ^ {(w + 2AK - vw x + (w + A) [{w + 2A) 2 - v 2 ] } , (4.6) 



where 



<t>x ' 



(4.7) 



w = (4-8) 



and A is an arbitrary constant. 
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• The singular manifold equations can be written as the system of PDEs: 

v t = (w x + wv) x (4.9) 

/ v 2 3 \ 

w t = \v x - — + -(w + 2X) 2 -2X(w + 2X)j . (4.10) 

Notice that the parameter A can be removed from these equations through the galilean 
transformation 

w -> (w + 2A), 
x — > x — 2Xt 

that transforms (4.9-10) into 

v t = {w x + wv) x , (4.11) 

m=(v x -^ + ?-w 2 ) . (4.12) 



2 2 , 

This system is equivalent to the Kaup or classical Boussinesq system. In fact it could be 
written as a single equation if we set 

w = p x . 

One can now remove v from (4.11-12) and the result is a sort of modified Boussinesq 
equation 

Ptt ~ Pxxxs - ^PxPxt + QplPxx - 2p t p xx = 0. (4.13) 

Unlike the AKNS case in which the Singular Manifold equations were also AKNS, for 
NLBq the Singular Manifold equations (4.11-12) are not the same system but rather they 
become the Kaup system (KS). As we will see in the next section this means that both 
systems NLBq and KS are related by a Miura transformation. 

4.2 Lax pairs 

The question of the linearization of the singular manifold equations is now a little bit more 
complicated that it was for AKNS. If we look for the dominant terms of (4.11-12) 

v ~ v x a , 
w ~ woX b - 
The result is: 

a = -l, b = -1, v = Xx, w = ±Xx, 



the ± sign confirms the well known fact that the Kaup system has two Painleve branches 
U, |§. In the previous section we have explained that for those types of systems it 
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is necessary the introduction of two Singular Manifolds. From that point of view, the 
truncation ansatz for v and w is now: 

»=^ + ^^fe=* + «-. (1.14) 

or 

2^f = v + w + 2A, (4.16) 

2^ = v-w-2X. (4.17) 

With this ansatz the expressions (4.5-6) for the truncated solutions (see Appendix D) 
can be linearized as: 

= 2M x ty+ x + M x i, + ) - (M t + M xx + 2AM X )^+, (4.18) 

= Vi + - ^+ + 2A^+ - 2M^+ (4.19) 

= 2M X (^" + M x ip~) + (i^ - + 2AM X )V>", (4.20) 

= + + 2AV^ + 2M x i>~. (4.21) 



To summarize, the existence of two branches in the Singular Manifold equations implies 
the existence of two classes of eigenfunctions of two different Lax pairs (4.18-19) and (4.20- 
21). This result reflects the fact that the equation is invariant under the discrete symmetry 

t — ► -t, 
M — ► —M 

that transforms the first of the Lax pairs into the second one. 



4.3 Darboux transformations 

The generation of the Darboux transformations can be done in a similar way as we did in 
Section 2. If we write the Lax pairs for the iterated solution M' 






= 2M' X ^' 2 + XX + M'J 2 + ) - (M' t + M' xx 




(4.22) 





= 4|-^ + 2A 2 ^+-2M>;+, 




(4.23) 





= 2M'M 2xx + M'Jf) + (Ml - M' xx 


+ 2X 2 M x )^ 2x , 


(4.24) 





= v47 + + 2A 2 4; + 2M>;-, 




(4.25) 
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where if) 2 + > ^2 cou ld be related to a singular manifold 2 in the form 

4>' 2 = ^' 2 -. (4.26) 

The consideration of (4.22-25) as coupled nonlinear equations allows us to write truncated 
expansions for M', t/j 2 + ,t/j 2 ~ , 2 . 

M' = M + %, (4.27) 
01 

4 + = V> 2 + + ^ (4-28) 
01 

1>2~ = i>2 + (4-29) 
01 

2 = 02 + I 1 , (4.30) 

01 

where 

= 2M X (^+ X + M^+) " (M t + M xx + 2\ 2 M x )il>+, (4.31) 

= ^ " 1>&x + 2A 2 ^+ - 2M^ 2 + , (4.32) 

= 2M X (^ 2xx + M x ip 2 ) + {M t - M xx + 2\ 2 M x )^ 2x , (4.33) 

= V2I + + 2A 2 fe + 2M^ 2 " , (4.34) 

02x = ^2^2" > (4-35) 

= 2M x (V'i a; + M x ^+) - (M t + M xx + 2A 1 M E )Y>+ , (4.36) 

= Y>5 - ^ + 2\ 1 *Pt x - 2M x ^t, (4.37) 

= 2M X (^ XX + + ( M t ~ M xx + 2AiM x )^, (4.38) 

= V'll + + 2Ai^- + 2M^r, (4-39) 

<f>ix = ipf^. (4.40) 
By substituting the truncation (4.27-30) in (4.22-26) it is possible to obtain 

e+ = -4jfn+, (4.4i) 

0- = -ip^n-, (4.42) 

A = -Q + tt-, (4.43) 

where 

n+ = (a^at) |f ( ^ 2+ * " (444) 

1 \ ^2 /„,,+ „,,+ _ „/, + „/,+ 



= ( J W ^ " ^ + ^iJ- (4-45) 
(4.27-30) together with (4.41-45) define Darboux transformations for NLBq. 
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4.4 Hirota's function 

The generation of a new iteration can be done by using 0' 2 as the Singular Manifold for 
M' in order to construct 

M" = M' + ^k = M + ^, (4.46) 

02 r 12 

where 

ri2 = 0' 2 0i (4-47) 
and by using (4.30) and (4.43) 

ri2 = 02^1 ~n + n-. (4.48) 

4.5 Solitonic Solutions 

The easiest nontrivial solutions can be obtained from the seminal solution 

M = a x. (4.49) 
For this solution, exponential solutions of (4.22-25) are 

ipf = exp{aj(x - ait)}, ip~ = exp \ -— ( x - — t ] > , (4.50) 



where Oj are related with the spectral parameter in the form 

Xi = ai + ^ (4.51) 

ai 

and the corresponding singular manifolds are 



au 



(ai + ^r), (4.52) 



where aij are arbitrary constants. (4.44-45) gives: 



./.-./,+ o- 01 

2<: 

and (4.48) is: 



o+ = * — n- = ^jty 2 ~ (4-53) 

a 2 ai - a a 2 ai - a 



Tl2 = 7^ ~T?2 V ("l + ^l" ^1 )(«2+V'2 f V'2 )-7 ai ° 2 ^ V^l ^2 -( 4 -54) 

(af - a )(a| - a ) (a 2 ai - a ) 2 
We can write the first (Fig. 3) and second iteration (Fig. 4) as: 

M' = a x + — , (4.55) 

01 

M" = a x + — , (4.56) 

Tl2 

where 

0i = -^-(l + Fi), (4.57) 
af - a 
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Tl2 



a\a 2 a\a2 



(a\ - a )(a| - a ) 
(Xi = exp 



{1 + F 1 + F 2 + A^F^} 



(4.58) 



a 

&i ^Oi ( i 



F 



a; 



exp 



«i H * - ^oi 



t4i2 = a 



/ a2 - ai 
laia 2 - a 



(4.59) 
(4.60) 



5 The Kaup system. Miura transformation to NLBq 



In the previous section, we have seen how the Kaup system (KS) arises as the Singular 
Manifold equation for NLBq. That suggests a Miura transformation between KS and 
NLBq [26], [22|. Let us write the Kaup system in the form 

u t = r] xx + 2uu x , 
Vt = Uxx + 2ui] x . 
Note that if we set u = p x , (5.1-2) can be expressed as follows 



Ptt ~ Pxxxx ~ ^PxPxt + QP x Pxx - 2<PtPx 



0. 



(5.1) 
(5.2) 

(5.3) 



This is precisely the Singular Manifold equation (4.13) for NLBq. 
If we use Painleve series for u and r] 



0=0 



3-a 



7] 



j=0 



the dominant terms yield: a = b = 1 and 

= ±Xx, 
r)o = Xx- 



(5.4) 
(5.5) 



The existence of two branches in the Painleve expansion suggests the following change of 
functions: 



u = m — m, 
j] = m + rh. 

With this change the addition and subtraction of (5.1) and (5.2) yields 

m t = m xx + 2(m - m)m x , 

rht = —iTT'xx + 2(m — rh)rh x . 



(5.6) 
(5.7) 

(5.8) 
(5.9) 
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5.1 Miura transformation 

From (5.8) we can obtain rh. By substituting it in (5.9), the result is: 

ml(mtt ~ ni xxxx ) = 4m 3 x m xx + 2m x (m t m xx - m xx m xxx ) - m xx {m 2 t - m xx ) (5.10) 

that is precisely the (4.3) NLBq. In a similar way m can be obtained from (5.9). Its 
substitution in (5.8) takes the form 

rh 2 x {rh a - rh xxxx ) = Am x rh xx + 2rh x (rh t rh xx - rh xx rh xxx ) - rh xx {rh\ - rh 2 xx ) (5.11) 

that is again NLBq. In consequence, the splitting (5.6-7) leads to the possibility of con- 
structing Soliton Solutions of KS by linear superposition of two solutions of NLBq. Actu- 
ally the inversion of (5.6-7) can be written as: 

2m x = u x - u 2 + d x l u t , (5.12) 
2rh x = —u x — u 2 + d~ x ut (5.13) 
that is a Miura transformation between KS and NLBq [26]. 



5.2 Backlund transformations 

Although two solutions m and rh of NLBq can be used to construct (by means of (5.6-7)) 
a solution of KS, these solutions are certainly not unrelated. Actually (5.8-9) establishes 
the correspondent relation between m and rh. This relationship can be written in the 
form: 

rht + rh xx 

m = m H , (5.14) 

2m x ' v ; 

m = m-\ (5.15) 

2m x 

which obviously is the Backlund transformation that relates the two solution of NLBq. 

5.3 Two Singular Manifolds 

The Singular Manifold approach derived in the previous section can be applied to m and 
rh. The Painleve expansion takes the form 

, , 4>X k , , <Px 4>X , K 1R s 

m =m+ — ^^>u =u+— —, (5.16) 



rh' = rh + ^^ri' = r!+-^ + -^. (5.17) 

(j) 9 <p 

Noneless the Backlund transformations (5.14-15) imply that <j) and <p are not unrelated. 
The substitution of the Painleve expansions (5.16-17) (or alternatively in (5.1-2)) gives 
rise to (see Appendix E) the coupling condition 

%i=^% + i4, (5.18) 
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where 

v — w 

A = ^- + u, (5.19) 

i = ^-«, (5.20) 

and 

A = A, (5.21) 

V + w — V + w 

u = \ + . (5.22) 

By using the definitions (4.16-17), the expressions (5.19-22) are: 

A = u + X + ^, (5.23) 

W 

A = -(u + \) + &, (5.24) 

U + A = $L-^2L. (5.25) 
5.4 Lax pair for KS 

The derivative of the coupling condition (5.18) with respect to x is: 

A At- A J> V-W-V + W tp~ 

A x = A{v- A- A) = = ^L_^, ( 5 _26) 

z tp ip 

A At A 1\ V + W - V - W V>r i>t /r o^N 

A x = A(y — A — A) = = ^~P (5-27) 

These expressions can be easily integrated as: 

A = a^, (5.28) 

ip 

tp + 

A = b%-. (5.29) 

tp + 

Combined with (5.23-24) these formulae yield 

ip~ = aip~ — (u + X)tp~, (5.30) 

^+ = aip + + (u + X)ip + . (5.31) 

The substitution of (5.25) in (5.30-31) leads to 

aip~ = (u x -rj x )ip~, (5.32) 
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aipt = -{u x + Vx)i> + , (5.33) 

where we have used i>^^ x = ~~ "^xwi>~-, i^t^x = —m x ip + ip~ (see Appendix D). 
Those expressions can be written as 

-(u + A) a 



U x Vx 







-iu x + rj. 



a (u + A) 




(5.34) 



(5.35) 



(5.34-35) are the spatial part of two components Lax pair for KS. The temporal part can 
be obtained from (4.17) and (4.19). 



t/>- 



4>~ 



^ o 1 iVxx - Uxx ~ (u - X)(rj x - u x )], Ux — x 
2a J 2 



V 

/ U x +T] x 

V 



Vx + Ux +U 2_ X 2 



a[u — A), 



2a [ ~" r " " " 



a{u — A) 



(u - X)(j] x + u x )] 



1>- 

4>- 



Vx^x +U 2_ X 2 



(5.36) 



.(5.37) 



6 Conclusions 



This paper has been dealing all along with the Painleve analysis in the version formulated 



by Weiss, Tabor and Carnevale |35]. Even though there is no rigorous proof so far available 
of the connection between the Painleve property and integrability the work hereby reviewed 
aims to contribute to a better understanding of the validity and usefulness of methods 
based on the Painleve property for studying Nonlinear Partial Differential Equations. 
With this idea in mind we would like to underscore some of the results that we have 
obtained in this paper by applying the Singular Manifold Method of Weiss [53] 



• In Section 2, we have applied the SMM to AKNS. This method have been proved 
to be quite useful to construct the Lax pair of AKNS. By applying the SMM to the 
Lax pair itself Darboux transformations and Hirota functions have been constructed 
algorithmically. The use of the SMM to construct solutions iteratively has been 
shown with the help of examples. 

• A similar procedure has been used in Section 3 to study NLBq. 

• The identification of Miura transformations and Backlund transformations by means 
of the SMM for equations with two Painleve branches appears as the main goal 
of Sections 3 and 5. There sinh-Gordon and Kaup systems are presented as the 
modified versions of AKNS and NLBq respectively. Two component-Lax pairs for 
both systems are obtained from the AKNS and NLBq Lax pairs as induced by the 
Miura map. 
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A Appendix 

By substitution of (2.3) in (2.1) we obtain a polynomial in ( — J whose coefficients are 

\<P J 

(we have used the code MAPLE V for the algebraic computer algebra): 
• Coefficient in ( — | 




AM y + 2v y + 8qM x + qv 2 + 2qv x = 0. 



(A.l) 




-4M xy - 2qM xx - 6vM y - \2qvM x - 8q x M x - 2q x v, 

2 7 1 3 , 

-2q x v - -vv x q - -qv xx - 3vv y - -v =0. 



•x 



(A.2) 




8vM xy + 4M xx (q x + qv) + AM y {v x + v 2 ) + 8M x (qv 2 + vq x + v y ) + v yxx 
+3vv xy + v xx (q x + qv) + 3v y (v x + v 2 ) + 3qv 2 v x + 3vv x q x + v 3 q x + qv 4 = 0. 



(A.3) 




= M, 



yxxx 



+ AM y M xx + 8M X M 



xy 



(A.4) 



(A. 4) means that M is a solution of AKNS and (A.l) can be used to obtain 



M y = -2qM x - -v y - -qv 2 - -qv x . 



(A.5) 



The substitution of (A.5) in (A.2) is 




that can be integrated as 



M x = 



v x v 2 X(t) 



(A.6) 



4 8 2 
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where A is a constant for the integration with respect to x. The substitution of (A. 6) in 
(A.5) is: 

M y = + \(t)q. (A.7) 
The cross derivatives of (A. 6) and (A.7) yield 
v xy ~ vv y = 4\q x , 

dt 

that are the Singular Manifold equations. 

B Appendix 

• Let us substitute (3.22) in (3.3). The result is: 

^ v + 2 u) + k(v-2u)-2^h = Q (B.l) 

<t> </> <t> <t> 

that compared with (3.24) reads 

A= V - + u, (B.2) 

A= V --u. (B.3) 

• Taking the derivative of (3.24) with respect to x. 

C Appendix 

The substitution of (4.4) in (4.2) leads to a polynomial in — whose coefficients are: 

3 



• Coefficient in , 

4M XX - w t + v xx + ww x + vv x = 0. (C.l) 
Coefficient in f — 



(C.2) 



N t - M xxx + 6vM xx - 2M t w - 6M% + M x (w 2 - v 2 - Av x ) 
+wi — v x + v v xx + v 2 v x + vww x — vwt = 0. 



Coefficient in 

vN t - vM xxx + 2M xx {v x + v 2 ) - 2M t (w x + wv) - QvM 2 
+M x (—v xx — v 3 + wt — 3vv x + ww x + w 2 v) = 0. 



(C.3) 
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If we set wt = ( v x — — + pj , equation (C.l) can be integrated in x as: 



1 



w 



p-v 2 -^- + 2X 2 (t) , 



4 V 2 

where A is an integration constant. By multiplying (C.3) for v and subtracting 

2w x (-2M t - vww x + 2wM x ) + (v x + 2M X )(AM XX - 2vv x ) = 0. 
The substitution of (C.4) in (C.5) yields 
1 



(C.2) 



-Avw x + 2(z - w)v x + (z + w) [ p - v 2 - ^- + 2A 2 (t) 



where we have equated 
p x = zw x . 

The compatibility M xt = M tx between (C.4) and (C.7) implies that 
z = 3w + 4A, 
3 

p= -w 2 + 4Xw + 2A 2 , 



dX 
~dt 



= 0, 



and therefore we obtain 

w t = (v x - V — + ^w 2 + 4Xw + 2A 2 

M x = - [(w + 2A) 2 - v 2 ] , 



Mt = ^{(w + 2X)v x - vw x + {w + X)[(w + 2A) 2 - v 2 }} . 

D Appendix 

To simplify the calculation let us define (see (4.16-17)) 



2a + = 2-f±- = v + w + 2X, 
2a~ = 2-^L = v - w - 2A, 



(C.4) 
(C.5) 
(C.6) 

(C.7) 

(C.8) 
(C.9) 



(CIO) 

(C.11) 
(C.12) 



(D.l) 
(D.2) 



or 



v = a + a 



(D.3) 
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w + A = a — a . 



(D.4) 



The substitution of (D.3-4) in (4.5-6) is 



M x = -a + 



a 



M t = -a + a- 2a+ - 2a~ - 2A + - 



or a 



(D.5) 
(D.6) 



In order to remove a from (D.5) and (D.6), we use (D.5) to set 



a. 



a 



+ 



(D.7) 



Its substitution in (D.6) provides: 



M t = M x 
= M 3 



.ot 



+ 



M xx M x , 
■— ^ + 2—^ + 2a + + 2^- - 2A 
M T a+ a+ 



M r . 



Mr 



ib + ib + 
+ 2M x ^ + 2^ -2A 



+ 



(D.8) 



that is (4.18). (4.20) can be obtained in the same form by removing a + between (D.5) 
and (D.6). 

• The temporal part of the Lax pair is obtained from the derivation of (D.l) with respect 
to t 



2at = 2^)= Vt + Wt . 
The use of (4.9-10) yields: 



(D.9) 



2a+ = 



v 2 3 

w x + wv + v x - — + -{w + 2A) 2 - 2X(w + 2A) 



that with the use of (D.3) and (D.4) finally leads to 
af = [a x + q +2 — 2Aa + — 4a + a~] ;r . 

Removing a~ with the aid of (D.7) 

a + = [ a + + a + 2 _ 2\a + - 4a + 2M x ] x . 

Finally we can substitute (D.l) and integrate out in x as 



(D.10) 



(D.ll) 



(D.12) 



(D.13) 



that is precisely (4.19). The expression (4.21) can be obtained by repeating the same 
process with (D.2). 
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E Appendix 

• The substitution of (5.16) in (5.1) gives 

— (-to + v + 2u) + i^(w + v-2u)-2- 

The comparison with (5.18) yields to 

v — w 

A = u + 



A = -u + 



2 ' 
v + w 



The derivative of (5.18) with respect to x provides just like in Appendix B 

A x = A(v - A - A), 

A x = A(v - A - A). 
The substitution of (5.16) in (5.2) gives us: 

^ • 2 2 

- \—w x — wv + v x + v + 2uv + 2i] x \ H — — [— w x — wv — v x — v + 2uv — 2i] x \ 

2 /? \ 2 



(E.l) 

(E.2) 
(E.3) 

(E.4) 
(E.5) 



4> 



+ 



+2- 



[w - v - 2u + 2A] + 



[w + v-2u- 2A] 



(E.6) 



v — v + 



0. 



The use of (E.4-5) and (5.18) yields to: 



— \-w x - wv + v x + v 2 + 2uv + 2??^ + 2A(v - v + A-A)} 

+^-[-w x - wv - v x - v 2 + 2uv - 2rj x + 2A(v - v + A - A)} = 0. 

Setting to both coefficients and using (E.2-5) 

2rh x = r] x -u x = -2A(v - A) = 0, 

2m x = r\ x + u x = — 2A(v — A) = 0. 
The comparison between (E.8-9) and (4.5) leads to 

v w + 2X 
A =2 ± —2—> 
A= v_ w + 2X 

2 2 



(E.7) 
(E.8) 

(E.9) 
(E.10) 



that compared with (E.2-3) means that A requires the plus sign and A the minus sign and 
A = A (E.ll) 
which means that 

W—V J V + W 1 r 

u = A-\ — = -A^ — = -[w-v + w + v + 2A]. (E.12) 
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figure 1 

One soliton solution M' x for k\ = 0, 2, k 2 = 0, 3 7 a = 0, 144 



figure 2 

Two soliton solution M" for k\ = 0, 2, fc 2 = 0, 3, a = 0, 144 



figure 3 

One soliton solution M' x for a\ = 1.2 a<i = 1.4, a = 3 



figure 4 

Two soliton solution M" for a\ = 1.2 0,2 = 1.4, a® — 3 



